Topological Phase Transition in the v — 2/3 Quantum Hall Effect 
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The double layer v — 2/3 fractional quantum Hall system is studied using the edge state formalism 
and finite-size diagonalization subject to periodic boundary conditions. Transitions between three 
different ground states are observed as the separation as well as the tunneling between the two layers 
is varied. Experimental consequences are discussed. 
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I. INTRODUCTION 

Advances in semiconductor fabrication have made it 
possible to produce multilayer two dimensional electron 
systems that allow exploration of the effects of interlayer 
as well as intralayer correlations. The possibilities raised 
by the introduction of an extra degree of freedom into 
the standard picture of the fractional quantum Hall ef- 
fect were first examined by Halperirill in the context of 
spin-uupolarized ground states, and later by Haldane and 
Rezayin who proposed applying Halperin's wavefunctions 
to the case where the electrons possess a double valued 
index indicating the quantum state of the electron in the 
third direction parallel to the magnetic field. These states 
haveqbeen analyzed in their spin as well as layer index 
-—Recent experiments with multilayer electron sys- 
: 'U seem to suggest the existence of incompressible 
states which belong to the universality class described by 
Halperin's wavefunctions. These systems allow the ma- 
nipulation of two parameters which greatly influence the 
character of the fluids, the distance between the effective 
layers of electrons, given by d and the tunneling between 
the two layers, denoted by the energy difference between 
the lowest two subbands as A sas . 

Motivated by the experimental observation of a tran- 
sition between two distinct,^ = 2/3 states in wide quan- 
tum wells by Suen, et.al.El and also the transition be- 
tween spin-polarized and spin-singlet v =r-2/3 states in 
tilted field experiments by Eisenstein, et.alfl, we have an- 
alyzed the double-layer model with interlayer tunneling 
at total filling factor v — 2/3 using the edge state for- 
malism and a finite-size numerical diagonalization study. 
In the double layer model electrons are localized on a 
pair of parallel planes, between which they can tunnel. 
Tunneling lowers the energy of the electrons in a sym- 
metric combination of states in each plane relative to the 
anti-symmetric state. If symmetric states are identified 
as pseudospin "up" states, and anti-symmetric states as 
pseudospin "down" , tunneling acts as a Zeeman coupling 
in pseudospin space. Therefore our investigation of the 
double layer model will allow a very general picture to 
emerge of two component systems where the pseudospin 
can refer either to real spin or to subband index. 



We have found a rich phase diagram with three distinct 
phases, separated by what we identify as first-order tran- 
sitions: (1) The particle-hole conjugate of the v = 1/3 
Laughlin state, where all the electrons are in the sym- 
metric state of the double layer, (2) The pseudo-spin 
singlet analog of the v = 2/3 spin-singlet state, where 
the electrons are divided equally between symmetric and 
antisymmetric states, and (3) a state with independent 
v = 1/3 Laughlin states in each layer. An interesting fea- 
ture is that states (1) and (2), while distinct, ,-have been 
identified as having the same topological orderQ, different 
from that of (3) . We carried out our investigation in the 
periodic (or toroidal) geometry, which is well adapted to 
exhibit these differences. For example, the long distance 
effective Chern-Simons (CS) theory of the Hall effect pre- 
dicts a three-fold ground state degeneracy for states (1) 
and (2) in the periodic geometry, but a nine-fold degen- 
eracy for state (3), consistent with our observation. An 
interesting consequence of the fact that state (3) has dif- 
ferent topological order from states (1) and (2) is that 
the boundary between state (3) and either state (1) or 
(2) necessarily supports a residual neutral gapless Lut- 
tinger liquidcl. The coexistence of two, and possibly three 
phases of quantum Hall states at v = 2/3 raises many in- 
teresting questions, both experimental and theoretical. 



II. EFFECTIVE THEORIES AND EDGE STATES 

Rather than work with the wavefunctions themselves, 
we can alternatively discuss the fractional quantum Hall 
states in terms of the long-wavelength effective field the- 
ories describing the incompressible fluids. These low- 
energy effective theories capture the long-distance corre- 
lations between the particles, determining such universal 
quantities as the conductance and the charge and statis- 
tics of the quasiparticles. The Lagrangian density for the 
effective theory can be written in the form 
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iftfaAjK^a*) +eA A (q,d jU a„)) , (1) 

where is an n-componcnt vector of Abelian CS gauge 
fields, Afj, is the electromagnetic gauge field, K is a non- 
singular integer coupling matrix, q is an integer vector; 
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(a, b) is the inner product. We note here that we use the 
notation q for the charge vector rather than t as used 
by Wen and ZeeQ. By integrating out the Chern-Simons 
fields a„ we can determine the Hall conductance of the 
above theory to be 



(2) 



The effective theory allows vortex defects of the Chern- 
Simons fields, with core energies determined by short- 
distance terms in the Hamiltonian which haven't been 
included. One can determine the charge and statistics of 
a particular composite of vortices by specifying an inte- 
ger valued vector 1 such that the composite in question 
is made up of vortices of type i. The charge of this 
composite is then given by 



(3) 



and the statistical phase acquired when two such com- 
posite vortices are interchanged is given by 



- = (1,K- 1 1). 
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Furthermore, it has been shownB that the degeneracy of 
a state described by the matrix K when it is defined on 
a two-dimensional closed surface of genus g is given by 
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(5) 



We may therefore classify an Abelian quantum Hall state 
by specifying an integer valued pair {K, q}0, thereby de- 
termining the long-distance properties of the fluid. It 
is important to note that distinct quantum Hall states 
are represented by equivalence classes of {K,q} pairs as 
the above properties are invariant under SL(k, Z) basis 
changes 
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where W is an integer matrix with |DctW| = 1. Two 
fractional quantum Hall states described by {K, q} pairs 
which are related by an SL(k, Z) transform belong to the 
same universality class and are considered topologically 
equivalent. 

The above discussion has been limited to Abelian frac- 
tional quantum Hall fluids on the plane. If we instead 
define our effective theory on a sphere, there is an extra 
term in the effective theory which describes the coupling 
of the liquid to the curvature of the sphere. This man- 
ifests itself in a new topological quantum number, the 
flux shift, which is not determined by the long distance 
effective theory described above. The shift S, defined by 
the relation = u N e — S, is a manifestation of the 
coupling of the orbital angular momentum properties of 
the state to the curvature of the space. The effective the- 
ories employed in the composite fermion approach on the 



plane do not distinguish between Landau levels and spin 
states. As the orbital angular momentum carried by by 
the cyclotron motion of electrons in the second Landau 
level is different from electrons in the lowest level, the 
shift S provides a way to distinguish between states that 
possess the same long distance properties but which have 
different spin symmetries. By utilizing the topological 
character of the coupling of orbital angular momentum 
to curvature we can classify the spin structure of a state 
whose analysis would otherwise lie outside the effective 
theory approach. 

We can write the Halperin-Laughlin wavefunction ap- 
propriate for multi-component systems in the planar ge- 
ometry in the form 



*K({^})=n^-^) KO ' (CTi,CT3) ii ex p(-^2 



i<3 



(7) 



where Gi is the pseudospin variable and Kf- = K el {ui,Uj) 
is the symmetric matrix encoding the electron correla- 
tions. If we specialize to the double layer system, we can 
identify the {mi,m,2,n) state with the effective theory 
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We shall be considering two effective theories for the 
system at v = 2/3. The first has been identified as rep- 
resenting both the pseudospin singlet state and the pseu- 
dospin polarized, particle-hole conjugate of a Laughlin 
v = 1/3 state. The effective theory is given by 
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The second effective theory that we will consider repre- 
sents two independent v — 1/3 Laughlin states 
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Both effective theories potentially represent states at 
v = 2/3. There is however a crucial difference between 
the two theories: they possess different ground state de- 
generacies on a non-trivial closed space, and therefore 
possess different topological order. In the remainder of 
this paper we will investigate the consequences of the 
fact that the these two theories have different topological 
order. 

Let us consider the edge state theory of the Abelian 
quantum Hall states. The edge between Abelian Hall 
fluids is a one-dimensional "Luttinger liquid" which can 
also be characterized by a {K, q} pair, which is the same 
pair of the bulk theory if the edge is between a Hall state 
and a non-conducting state. Generally, we may define a 
set of fields <fii(x) living on the one-dimensional compact 
edge of an incompressible Hall sample such that they 
obey the equal time commutation relations 
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[<fii(x), 4>j{x')} = iir(K i:j sgn(x - x') + L i} ) (11) 

where Lij = sgn(i — j)(Kij + qiqj) is a Klein factor. The 
action for the translationally invariant edge fields is given 

by 



/ 



where 



and 



S = dt(S - H ) 



S = ^2 K ij 1 j> dx 



H ° = ~^ Yl Vi i j> dx d x4>id x (i 



(12) 



(13) 
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where Vy encodes the non-universal interactions which 
determine, among other things, the velocities of the var- 
ious edge modes. For stability, we require the matrix V 
to be positive definite. Let us define 



and 



q(m) = (m,q) 



(15) 
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where m is an integer valued vector. We can define a set 
of composite local fields 

* m = e - ^" (17) 

which obey, for x =/= x' 

^m(x)^ m '(x') = r}m, m >V m >(x')V m (x) (18) 

where 

Vm, m > = (-1)1(™M™'K (19) 

The winding number operator, defined as 



Ni = — (f> dx d x <f>i(x) 



(20) 



is constrained to be integral by imposing periodic bound- 
ary conditions on the fields 

V m (x + L) = y m (x). (21) 

The fields are characterized by the integer quadratic form 

K(m) = (m, Km) = m^rn,- (22) 



which has the constraint 



(23) 



If K(m) is odd, the field corresponding to m is is 
fcrmionic, and if K(m) is even, it is bosonic. We can 
define the charge density of the edge as 



which has the commutation relation 

[p(x),p(x')]=iha H 6'(x-x') 

where 
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(24) 



(25) 



(26) 



The physical interpretation of this quantity is that it rep- 
resents the change in the quantized Hall conductivity in 
going from one side of the edge to the other. We can also 
define the total charge operator 



ij 

which obeys the commutation relation 
[Q,V m (x)] = eq(m)^ m (x) 



(27) 



(28) 



The edge state Hamiltonian Ho describes n linearly 
independent oscillator modes propagating with velocities 
v\ which can be determined from the generalized real 
symmetric eigenvalue equation 



Vu A = v x K 1 u x 



(29) 



where the velocities are real and the eigenmodes inde- 
pendent as V is positive definite. We can now state the 
condition on the matrix K which ensures a set of SU (2) 
generating operators within the theory: If we can identify 
an integral vector m such that 



K (m) = 2s m q(m) = 



(30) 



where s m = ±1 then we may identify an SU(2) algebra 
associated with the edge field * m . We can define a triplet 
of non-Abelian densities 

a x (x) = i(* m (a;) + *_ m (z)) 

^) = TSf(*m(*)-*-m(aO) 

a z (x) = \p m {x) = -^d x <t> m (x) (31) 
which obey a level- 1 SU{2) Kac-Moody algebra 



[a°(x),a b (x')} = t -^-6>(x-x>) 

+ ie abc a c {x)S(x - x'). 
We may identify the SU{2) algebra generators 



S a = 



j> dx a^(x) 



(32) 



(33) 
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which obey the commutation relations 

[S a ,S h ]=ie abc S c . (34) 

The auxiliary constraint q(m) = is to insure charge 
neutrality.r,T4iis hidden SU (2) symmetry has been noted 
previouslytHlil. We note that larger symmetries (specifi- 
cally SU (N) for an N-component system) may be iden- 
tified using a similar analysis. 

Let us focus specifically on the double component sys- 
tem, where the vector components denote different pseu- 
dospin components. We shall work in the symmetric ba- 
sis where % = 1 for i = 1,2. Further, the diagonal ele- 
ments Ka must be odd for fermi statistics. We therefore 
find that the only effective theory which may generate an 
SU{2) edge symmetry is 

K=( m±1 ™) (35) 
V m m ± 1 J 

with m even. Diagonalizing the above edge state sys- 
tem ( |35"| ) we obtain the Hamiltonian 

H = |~ / dx i v n{dx<l>n(x)) 2 + V c (d x (j) c (x)) 2 

+2v int d x (f> c (x)d x <f> n (x)) (36) 

where </>„ = 4>± — (f>2 is a neutral mode and and <p c — 
4>\ + 02 is a charged mode. The two modes move with 
velocities that depend on the external potential, with in- 
termode interactions Vi nt mixing these two modes. It is 
the neutral edge mode that generates the SU (2) algebra 
where the field 

*m = e~ iWl ~* a) (37) 

corresponds to the neutral physical operation of tunnel- 
ing one electron from one pscudospin to the other. The 
effective theory defined by (|3^) possesses a hidden SU (2) 
symmetry of the bulk quantum state, in this case in- 
variance under pseudospin rotations. We may therefore 
state that a bulk quantum Hall state may possess an 
SU(2) symmetry if and only if the effective theory cor- 
responding to this quantum Hall state has a non-trivial 
solution m to the equation K(m) = ±2. As the bulk 
ground state is non-degenerate, apart from topological 
degeneracies which are not associated with the SU(2) 
symmetry, our states constructed in this way are SU (2) 
singlets. Whether or not a particular state is realized as 
the ground state depends, of course, on the underlying 
bulk Hamiltonian. It is important to realize that this 
symmetry is a consequence of the electron correlations 
rather than the underlying Hamiltonian. This algebra is 
realized in the bulk state independently of whether or not 
the Hamiltonian is strictly invariant under the symmetry. 

We must note that in our construction we have as- 
sumed that the matrix K was non-singular. The singular 
case of Kij = m with m odd corresponds to a v = l/m 
state which is also invariant under pseudospin rotations, 



but which is fully polarized with S — N e /2. This system 
has been studied pr.6yiously and is found to possess many 
interesting featurese3. 

The edge state formalism is also capable of addressing 
the question of what happens at the edge between two 
Abelian Hall states, which is physically realized when 
there is phase coexistence in a first-order transition be- 
tween the two states. We can form the edge state theory 
corresponding to the {K, q} effective by forming the di- 
rect sum K = Ki © — K 2 and q = qi © q 2 . We will 
restrict our discussion to a clean edge between two dou- 
ble component Hall states at the same filling fraction, at 
least one of which possesses an SU(2) symmetry of the 
type discussed above. To address the stability of the edge 
of this system, we must consider the general non-winding 
number conserving tunneling perturbation 

T = J dx[t(x)^ m (x)+c.c.]. (38) 

The operator must be bosonic and charge conserving with 
q(m) — and K(m) even. The scaling dimension of this 
operator can be determined from the two-point function 
to be 2 — A(m), where A(m) satisfies the inequality 

A(m) > ±\K(m)\. (39) 

We can construct a representation of the fields such that 
both Kij and the non-universal Vij are diagonal 

{UKU T \, = cnSij (40) 

(C/y(7 T )y = ViSij (41) 

where Ui = ±1. Since stability requires V{ > 0, the di- 
rection of propagation of each mode is determined by 
a - , . Using the transformation matrix U we can determine 
A(m) 

5^7fc(m)=2A(m) (42) 

i 

where 

1.H = E^ 1[/ hS™^ (43) 

The factors 77(m) obey the sum rule 

^^(m)^ = K(m). (44) 

i 

One can see that if the system is maximally chiral, then 
the sum rule explicitly gives us the scaling dimension of 
the operator. If it is not maximally chiral, it only gives 
us a lower bound. Therefore, from a scaling perspective 
a tunneling operator with K(m) = is potentially rel- 
evant if the scaling dimension A(m) < 2. In general, 
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the tunneling perturbation will be prevented from being 
relevant by the complex tunneling parameter t{x) where 



t{x) = \t(x)\exp[ia{x)] 



(45) 



in the clean case. In order for the perturbation to be 
relevant, the phase factor must satisfy 



d x a(x) = (0\d x (j> m \0). 



(46) 



We will assume that such a "phase locking" is generi- 
cally possible for the edges in question. If the K(m) = 
perturbation is relevant, the modes involved in <I> m be- 
come massive and are removed from the low-energy the- 
ory. The only other perturbations that can be poten- 
tially relevant have |if(m)| = 2 but are not mass gen- 
erating (they do form 'hidden' SU (2) symmetries as dis- 
cussed previously). Therefore, the condition for a po- 
tentially mass generating perturbation is that we must 
identify a non-trivial integer valued vector m such that 
K(m) — q(m) — 0. We shall take the point of view 
in this paper that if a mass-generating perturbation is 
potentially relevant, the system will relax so that the 
instability generally occurs. This is a conjecture based 
upon the experimental observation that the only stable 
Abelian Hall states are those that do not permit such 
mass-generating instabilities^. If an mass-generating in- 
stability can occur, it appears to do so. 

We can apply this analysis to the edge theory at hand. 
First, let us consider the edge between two states with 
the same topological structure. The edge state theory 
will be based on the {K,q} pair 
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(47) 



(48) 



We note that the above edge has <jh = as is appropriate 
for an edge between two Hall states with the same filling 
fraction. In this case, we can identify the two operators 
of interest: the operator ^ m with m = (1,0,-1,0) and 
the operator with m = (0, 1, 0, —1). Both these operators 
are K(m) = operators which are potentially relevant 
and mass-generating. Interactions will tend to reduce 
the scaling dimension of the two operators, but if they 
remain relevant we expect that they will cause a mass gap 
to form (we note that generally diagonal elements of the 
interaction matrix reduce the scaling dimension of the 
operators, while off-diagonal elements tend to increase 
it, upto a maximum of 2 — ^\K (m)|). It is straightfor- 
ward to see that any edge between states with the same 



K matrix will have two such tunneling operators with 
K (m) = q(m) = which potentially cause the modes 
to pair up and form a gap, leaving no residual gapless 
modes in the low-energy theory. Let us consider a single 
quantum Hall state where we arbitrarily choose there to 
be an edge in the bulk of the state. The edge state will 
have the same K matrix structure as the edge between 
any two states that have the same topological structure. 
Such a fictitious edge in the bulk cannot support gap- 
less charged excitations, as is consistent with our above 
analysis. Generically, we then expect that at the edge 
between two states with the same topological structure 
there should be no residual gapless modes. 

We can also consider the edge between two states with 
different topological structure at the same filling fraction, 
as is appropriate at the edge between the pseudospin sin- 
glet state or the particle-hole conjugate state and the 
state with independent v = 1/3 Laughlin states in each 
layer. Generically there is only one K(m) = tunneling 
operator at the edge between double component states of 
different topological order. It isn't possible for two sets 
of edges to pair up and form a gap, due to the mismatch 
in topological structure. There will always be a pair of 
residual neutral gapless edge modes left over in the low- 
energy theory, with a gap for making charged excitations. 
As an example we shall consider the edge between two 
states at v — 2/3 possessing different topological order. 
This will be based on the {K, q} pair 



K = 



1 


2 
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1 
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1 -1 




















(49) 



The operator ^f m with m = {1, 1, — 1, — 1} has K(m) — 
and g(m) = 0, thereby allowing two charged modes 
to pair off and form a gap, leaving two neutral gapless 
modes in the low-energy effective theory. The mismatch 
between two distinct quantum Hall states at the same 
filling fraction with different topological order implies the 
existence of residual neutral gapless excitations at the 
boundary. 

In considering the edge between two distinct quantum 
Hall states at the same filling fraction, we expect two pos- 
sible scenarios. If the two states possess the same topo- 
logical structure, such as the pseudospin singlet state and 
the particle-hole conjugate state, their respective edges 
will pair up and form gaps, leaving behind no residual 
gapless states. If the two states have different topological 
structure, only one set of modes will split off and form 
a gap. Two neutral, gapless modes will remain in the 
low-energy theory, a residual side effect of the mismatch 
in topological order. 
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III. PHASE DIAGRAM FROM FINITE SIZE 
STUDIES 

The finite size studies that we report were carried out 
using periodic boundary conditions (see Appendix) and 
delta function wavefunctions to represent the two layers 
in the double layer model, with the electrons confined to 
the lowest Landau level as is standard practice. Our stud- 
ies were performed at filling fraction v — N e /N$ = 2/3 
with six electrons and nine flux quanta. In the following, 
length will be measured in units of I, the magnetic length 
^jTi/eB, and energy in units of e 2 /Airel. The system of 
electrons interacting via the Coulomb interaction is ex- 
actly diagonalized numerically, with the spectrum as a 
function of k (see Appendix) providing the fundamental 
information on the system. We have diagonalized the 
system with the Hamiltonian given by 



H = - 



^■sas 



\ ^{^,jaj\V\ka k ,lai)c\^ aj c kak ci ai (50) 



where 



i,j,k,l 



V(r) 



4ne (r 2 +rf 2 (l-«5 Ql ^)) 1 /2 



(51) 



and Ui is the layer index which denotes in which of the 
two layers the electron resides. 

All calculations were performed using square boundary, 
conditions where 9 = n/2 and |Li| = |La|. It is knownES 
that incompressible states are remarkably insensitive to 
the particular boundary conditions chosen, as long as 
the shortest length scale of the geometry is larger than 
average interparticle spacing. While the exact details of 
the excitation spectrum in the system under investigation 
will depend on our choice of conditions, the qualitative 
conclusions concerning the incompressible ground states 
should not. 

In the following we shall use the term pseudospin to 
refer to the subband layer index, with the the up-spin cor- 
responding to the symmetric combination of layer states 
and the down-spin the antisymmetric. Tunneling acts as 
a Zeeman term in the Hamiltonian which tends aligns 
the electron pseudospin in the up state, or the symmet- 
ric combination of layer states. We will vary both the 
tunneling, denoted in the Hamiltonian by A sas , and the 
distance between the double layer planes, denoted by d/l, 
investigating both the ground state and the dependence 
of the excited states on these parameters. 



A. Spin-Singlet State 

It was realized some time ago0 that when the electron 
correlations are of the same order as the Zeeman energy 



associated with the spin states, it is important to consider 
the spin degrees of freedom in constructing the ground, 
state. Previous numerical and experimental studies&ta 
on the v = 2/3 system reveal that the ground state of 
the non spin-polarized system, in the limit of vanishing 
Zeeman energy, is plausibly a spin-singlet. If we consider 
a phase diagram where we vary both the interlayer sep- 
aration d/l and the interlayer tunneling A sas , along the 
line d = the Hamiltonian is invariant under pseudospin 
rotations, allowing a direct mapping between pseudospin 
and electron spin in the presence of Zeeman term. We 
know from our previous edge analysis ( |33| ) that the K 
matrix of the effective theory of a spin-singlet ground 
state must be of the form 



K 



777 ± 1 777 
777 777 ± 1 



(52) 



where m is even. Jainli3 has proposed a wavefunction to 
describe the spin-singlet state at v = 2/3 based upon the 
effective theory where {K, q} is given by 



K 













i 


q = 





(53) 



We cannot strictly interpret this K matrix in terms of 
the Halperin-Laughlin wavefunction, as thermodynamic 
stability requires the K matrix for such a wavefunction 
to be positive definite, so we must interpret the wave- 
function in the planar geometry to be 



'3 i<3 

II exp ( i f sgn ( ' i ~~ a i 

i<3 ^ 



l<3 

x 



(54) 



where we have denoted d Zi = . This wavefunction is a 
spin-singlet at the correct filling fraction with the correct 
topological properties. 

In Fig. |lj we calculate the overlap of the exact ground 
state with Coulomb interactions with Jain's proposed 
spin-singlet wavefunction as we vary the Hamiltonian pa- 
rameters A sas and d/l (note that A sas and tunneling are 
used interchangably) . Several points should be noted. 
First, the system is very sensitive to the effects of tunnel- 
ing, or equivalently, a magnetic field in pseudospin space. 
The overlap with the spin-singlet state rapidly falls off 
with the introduction of even a slight amount of tunnel- 
ing. This will tend to make the experimental observation 
of the spin-singlet state in multilayer samples extremely 
difficult. Second, the system is reasonably robust against 
a separation between the two layers, falling to zero at 
d rs 1.1/. Thus we find significant overlap between the 
exact ground state and the proposed spin-singlet even in 
regions where the Hamiltonian no longer commutes with 



G 



the pseudospin algebra. The spin-singlet character of the 
state is a manifestation of the electron correlations rather 
than the underlying Hamiltonian. 

While we expect that the numerical data gives us qual- 
itative data on the transitions discussed, it should be 
noted that finite size effects will influence the exact po- 
sitioning of the transitions in the d/l — A sas plane in re- 
lation to the thermodynamic limit, and should therefore 
be taken with caution. 



B. Particle Hole Conjugate of v — 1/3 Laughlin State 



The particle-hole conjugate of the standard Laughlin 
v = 1/3 gives us an incompressible liquid at filling frac- 
tion v — 2/3 which is denoted in the effective theor y b y 
the same K matrix as that of the spin-singlet state (53). 
The matrix nature of the effective theory does not re- 
flect the correlations between electrons of opposite pseu- 
dospin, as the particle-hole conjugate state is pseudospin 
polarized, but rather reflects its composite nature. While 
to our knowledge no simple wavefunction has been pro- 
posed for the exact particle-hole conjugate state, a trial 
hierarchy wavefunction which very effectively captures 
the electron correlations has been developedta. We may 
write the hierarchy v = 2/3 state as 



i' i'<j' 
i<j 

x JJexpf-^l^n 

i' ^ ' 

X Il eX p(-^2l Z 



* \2 



(55) 



where i' = 1, Nh and i = 1, N e where the number of 
holes Nh — N e /2 and N e is the number of electrons. This 
hierarchical state is represented by the effective theory 
given by 



1 1 

1 -2 



K= | , „ | q 



(56) 



which is equivalent to (^3() upto a similarity transforma- 
tion. This effective theory also possesses an SU(2) in- 
variance. In this case, the two pseudospin variables refer 
to the lowest two Landau levels in the composite fermion 
construction. Within the composite fermion approach, to 
obtain the single layer particle-hole conjugate state one 
starts with a spin polarized v = 2 state, a gapped sys- 
tem, and adiabatically attaches two flux quanta to each 
electron opposite to the direction of the magnetic field 
that generated the v — 2 state. As the addition of two 
flux quanta doesn't affect the statistics, our composites 
arc still fcrmions. In a mean- field sense, we start with a 



pseudospin polarized state at v = 2 and decrease the B 
field by two flux tubes per electron, arriving at a pseu- 
dospin polarized state at v = 2/3, which we identify as 
the particle-hole conjugate state. We still have resid- 
ual gauge fluctuations associated with the added flux, 
but they shouldn't qualitatively change the physics as 
we started with a gapped system. Therefore the v = 2/3 
polarized state can be identified with a polarized v = 2 
integer quantum Hall state. For the spin-singlet state, 
out starting point is a spin-unpolarized v = 2 state with 
the first Landau level being filled for both the up and 
down spins. We then perform the same flux addition 
process as we did for the polarized state to arrive at a 
spin-singlet v = 2/3 fermion state. Within this approach, 
which state will be realized depends on the ratio of the 
effective cyclotron energy to the effective Zeeman energy 
of the composite fermions. 

In Fig. H we calculate the overlap of the particle-hole 
conjugate state with the exact ground state as a func- 
tion of A sas and d/l (A sas and tunneling being used in- 
terchangably as in Fig. [Ej). The effect of the tunneling 
in the spin analogy can be seen to be a turning on of a 
Zeeman energy term in the z direction. As this energy is 
increased, eventually all the spins will align themselves 
along the z direction, resulting in a spin-polarized Laugh- 
lin state with all the electrons occupying a symmetric 
combination of the layer indices. As the Zeeman term is 
turned on, the system abruptly reaches a transition point 
where the the overlap of the ground state with the spin- 
singlet state falls to zero (Fig. jij), while the overlap with 
the Laughlin 2/3 particle-hole conjugate state jumps to 
close to unity. 

As the spin-singlet state and the particle-hole conju- 
gate state are described by the same pair {K, q}, they are 
described by the same effective theory on the plane and 
possess the same edge state structure. While the effec- 
tive theory approach captures many of the long-distance 
properties of the fluid, it doesn't classify the spin of the 
state, which is determined by the energy associated with 
the spin degree of freedom. As we turn on the tunneling 
A sas , we find a transition from a region where it is ener- 
getically favorable to put the composite fermions in the 
first two spin states in the lowest Landau level to a region 
where it is favorable to put them in the first two Landau 
levels. In both cases we have an SU (2) symmetry, in the 
one case between pseudospin states and in the other be- 
tween Landau levels. It is important to notice that within 
the composite fermion approach the effective cyclotron 
energy is heavily renormalized, as small variation in the 
effective field essentially makes it energetically favorable 
to place the electrons in two pseudospin-polarized Lan- 
dau levels, costing cyclotron energy but saving on Zee- 
man energy. The composite fermion process of attaching 
flux has the effect of enhancing the ratio of the effective 
energy associated with a pseudospin flip to the effective 
cyclotron energy. 
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C. (3,3,0) Double Layer State 



IV. FINITE SIZE STUDY OF TRANSITIONS 



At non-zero d/l, as the Hamiltonian no longer com- 
mutes with psuedospin rotations the ground state need 
not be an eigenstate of pseudospin nor need it possess 
an SU (2) symmetry as the particle-hole conjugate state 
does. A proposed effective theory for the ground state at 
large layer separation d/l is given by 











! 


1) 


q = 


(!) 



which has a direct interpretation as a (3,3,0) Halperin- 
Laughlin multi-component wavefunction representing in- 
dependent v = 1/3 Laughlin states in each layer. In 
Fig. H we calculate the overlap of the exact ground 
state with the (3, 3, 0) state as a function of A sas and d/l 
(note again that tunneling and A sas are being used in- 
terchangably) . The transition between the (3, 3, 0) state 
and the particle- hole conjugate state (Fig. 0) represents 
a competition between minimizing the coulomb energy 
between the electrons and the tunneling energy cost of 
localizing the electrons on independent planes. For any 
given tunneling A sas it is possible to find a large enough 
d/l such that it is energetically favorable to form the 
independent v — 1/3 states in each layer, costing in tun- 
neling energy but gaining even more in minimizing the 
electrostatic ground state energy. Equivalently, given any 
interlayer distance we can find a tunneling energy which 
will outweigh the favorable interaction energy associated 
with having the electrons as far apart as possible. 

D. Phase Diagram 



Whenever a system exhibits different macroscopic 
phases as a function of system parameters, it is natural to 
ask questions about the transitions between such phases. 
While finite size studies are unable to address questions 
about the thermodynamic features of such transitions, 
they are able to shed light on qualitative changes in 
structure that a bulk liquid may undergo in going from 
one phase to another. In order to investigate the tran- 
sitions involved in our double layer v = 2/3 system, we 
have exactly diagonalized our finite system using periodic 
boundary conditions as a function of the system param- 
eters A sas and d/l. We can follow the development of 
individual energy levels of the system as these parame- 
ters are varied by identifying the quantum numbers of 
the state in question such as parity under reflection and 
the translational quantum number k and noting that un- 
der adiabatic perturbation the energy levels should be 
smooth and connected in our finite system. In the fol- 
lowing figures we plot the energy levels of the system as a 
function of either A sas and d/l, at each step subtracting 
off the average energy of the system in order eliminate 
background energies of the system. The energy levels 
which are relevant to each transition are connected for 
emphasis during their development. 

When using periodic boundary conditions it is impor- 
tant to note that each eigenstate has a generic degeneracy 
associated with center-of-mass translations given by q if 
the filling fraction is v = p/q. The formalism used to 
classify states (see Appendix) extracts this degeneracy 
explicitly. Since we are always working at fixed filling 
fraction, this degeneracy will be unimportant. 



We can identify three separate phases for the v = 2/3 
system as a function of tunneling A sas and the layer sep- 
aration d/l. We calculate the phase diagram for the 
v = 2/3 system in Fig. Q from the overlap data by 
matching a system with a particular phase if the over- 
lap of the ground state with the characteristic state is 
greater than some cutoff, chosen to be .75 in Fig. |Q. 
One important feature to note is the existence of a triple 
point in our numerical studies at d ~ l.ll and A sas w .01 
where all three phases will be in coexistence. It is diffi- 
cult to extract the experimental parameters where such 
a triple point might occur from our finite size numerical 
data, as such quantitative information will be sensitive 
to finite size effects such as geometry and particle num- 
ber. We conjecture that the triple point where all three 
discussed phases will be in coexistence will persist in the 
thermodynamic limit, in the neighborhood of the physi- 
cal parameters suggested by our studies. The observation 
of such a triple point remains an interesting experimental 
possibility. 



A. Variation of d/l at A sas = 

In Fig. H|, the variation of the energy levels as a func- 
tion of d/l for A sas = is shown, with the average energy 
at each d/l subtracted off. The transition being witnessed 
is from the spin-singlet phase to the (3, 3, 0) phase as d/l 
is increased. At small values of d/l the ground state 
is well represented by Jain's spin-singlet state, with a 
well-defined energy gap to all excitations indicative of an 
incompressible phase. 

At large values of d/l the ground state is given by 
a three- fold degenerate multiplet of states, each being 
a k = eigenstate which in the thermodynamic limit 
becomes rotationally invariant. There is again a well- 
defined energy gap between the ground state multiplet 
and all excited states. There is a slight splitting of the 
degeneracy due to finite size effects but this feature will 
disappear as the size of the system is increased. 

The key feature to note is the change in ground state 
degeneracy in going from the spin-singlet state to the 
(3,3,0) state. Such a change is to be expected from 
our effective theory considerations. As noted previously 



8 



the degeneracy of a state given by the effective theory 
with the pair {K,q} on a torus is given by |Detif|. 
The spin-singlet state must then have an overall degen- 
eracy of three while the (3,3,0) state has a degeneracy 
of nine. As mentioned above our formalism extracts a 
three-fold center-of-mass degeneracy generic to states at 
v = 2/3, leaving us with a residual three- fold degeneracy 
for the (3, 3, 0) state and a non-degenerate spin-singlet 
state, consistent with our numerical data. 

Another point to be noted is that one of the energy lev- 
els coming down in the (3, 3, 0) state triplet has the same 
symmetry as the spin-singlet state, resulting in an energy 
level repulsion as they cross. As the transition point ap- 
proaches, the triplet of states comes down, crossing with 
the spin-singlet state. 

It is interesting to note that even in considering bulk 
transitions, as we are in our finite size studies, the sys- 
tems display a residual side effect from the topological 
mismatch between the two states in coexistence. In our 
edge state analysis we found that at the edge between 
the spin-singlet and (3, 3, 0) state there will be residual 
neutral gapless modes due to the difference in topologi- 
cal structure, while in our bulk finite size studies we find 
a change in ground state degeneracy. Both features are 
generic and stable against perturbation, indicative of the 
true topological character of the incompressible fluids. 

B. Variation of A sas at d/l = 

In Fig. B we investigate the energy level structure 
by fixing d/l = and varying A sas . We are witnessing 
the transition from the spin-singlet ground state to the 
particle-hole conjugate state as a function of the Zeeman 
energy associated with the pseudospin degree of freedom. 
The system quickly undergoes a transition from the spin- 
singlet state to the particle-hole conjugate state indicated 
by a simple level crossing. The spin-singlet state is ex- 
tremely sensitive to the effects of the effective magnetic 
field, rapidly finding it energetically favorable to place 
the composite fermions in the second Landau level rather 
than the first spin-reversed Landau level. 

As the two phases possess the same topological struc- 
ture they also possess the ground state degeneracy. They 
do however have different pseudospin symmetries, allow- 
ing the energy levels to cross without repulsion. Such an 
energy level crossing will become a first-order transition 
in the thermodynamic limit. 

C. Variation of A sas at d = 2.0Z 

In Fig. 0| we investigate the energy structure at d = 
2.01 as we vary A sas . We are seeing the transition from 
the (3, 3, 0) state to the particle-hole conjugate state as 
we turn up the tunneling. The transition is qualitatively 
the same as the spin-singlet to (3,3,0) transition, with 



a degeneracy transition due to the topological mismatch 
between the two states. Again we see a level crossing 
driven by variation in the sample parameters, where two 
states involved in the crossing mix and cause energy level 
repulsion. This transition is particularly relevant as ex- 
perimental evidence supporting such a phase transition 
already existsEl. The experiments of Suen, et.al. were 
performed using a single wide quantum well geometry, 
making quantitative comparison with our idealized dou- 
ble layer calculations difficultLj. 

V. NATURE OF THE TRANSITIONS 

We believe that the transitions seen in our finite size 
studies represent first-order phase transitions in the ther- 
modynamic limit at low temperatures in a clean system. 
Traditionally, one uses broken symmetries and their as- 
sociated order parameters to classify and organize many- 
body systems, often allowing the construction of low en- 
ergy effective theories based on these order parameters 
which capture the essence of the correlated states as well 
as predicting effects that are not accessible in a micro- 
scopic approach. The effective theories employed in the 
description of the fractional quantum Hall effect are not 
based on an order parameter derived from a broken sym- 
metry. Rather, these effective theories embody a new 
type of order called topological order which is stable 
against perturbation and manifests itself in such prop- 
erties as the ground state degeneracy when the system is 
defined on a topologically non-trivial closed space. We 
therefore do not consider the fractional quantum Hall 
states to be broken symmetry states. 

As the fractional quantum Hall states are not broken 
symmetry states, they cannot undergo the usual sec- 
ond order phase transition where the order parameter 
smoothly goes to zero in the vicinity of a transition point, 
and is zero on the other side of the transition. It would 
appear that due to the fact that the order embodied by 
the fractional quantum Hall states is discrete and topo- 
logical in origin that it would be impossible to smoothly 
go from one state to another with different topological or- 
der as a function of system parameters in a clean system. 
Let us consider a point in the d/l — A sas plane where 
either the v = 2/3 spin-singlet state or the particle- hole 
conjugate state is in coexistence with the (3,3,0) state. 
From our edge state analysis we know that when two 
phases having different topological order are in coexis- 
tence there will be at least one pair of residual neutral 
gapless modes at the boundary between phases. The fi- 
nite energy of these modes will localize the boundary 
between the two phases, forcing the transition to occur 
by nucleation, resulting in a first-order phase transition. 
If two states have different topological order, there is no 
way to smoothly go from one state to another due to the 
discrete nature of topological order. 

In the transition between two states that possess the 
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same topological order there is no such mismatch. If 
we consider the interface between two such states, such 
as the spin-singlet state and the particle-hole conjugate 
state, we expect that the two pairs of edge modes will pair 
up and form a gap, leaving no residual gapless modes. In 
a clean system the two states have different pseudospin 
symmetries, as the particle-hole conjugate state is spin- 
polarized while the spin-singlet state is a true pseudospin 
singlet. This difference is reflected in the flux- number 
shift of the two states on the sphere: 5 = for the 
particle-hole conjugate state and S = 1 for the spin- 
singlet. The transition can then go by a simple level 
crossing, resulting in a first-order transition in the ther- 
modynamic limit. 

There exists the possibility that there might be some 
intermediate state that exists between the two principal 
states undergoing the transition. This case is really a two 
step process rather than a direct transition: Hall state to 
intermediate state, then intermediate state to Hall state. 
In principal the intermediate state could be incompress- 
ible, but we can then apply the same arguments used 
above to show that the two principal states undergoing 
the transition cannot be smoothly connected. While this 
scenario is an experimental possibility, in the v = 2/3 
system no evidence exists for this type of intermediate 
transition. As such, we conjecture that all the transi- 
tions involved in the v — 2/3 system will be first-order. 



VI. CONCLUSIONS 

In conclusion, we have examined the structure of the 
phase diagram of the v = 2/3 double layer electron sys- 
tem as a function of d/l, the distance between the layers, 
and A SQS , the tunneling parameter for the system. A 
phase diagram consisting of three different phases, each 
belonging to a different, distinct universality class, was 
calculated. A triple point is conjectured where all three 
phases are stable. A gapless, neutral Luttinger liquid 
structure is predicted at the interface between either the 
spin-singlet state or the particle-hole conjugate state and 
the (3, 3, 0) state. At the interface between the spin- 
singlet state and the particle-hole conjugate state no 
residual gapless modes are expected. It is conjectured 
that there should be first-order transitions between all 
three phases, indicated by distinct level crossings in the 
finite system energy levels. 



APPENDIX: PERIODIC BOUNDARY 
CONDITIONS 

We wish to impose generalized boundary conditions 
by requiring that all physical quantities be invariant un- 
der translation of any particle by the set of translations 
L mn = toLi + nL 2 such that 



|Li x L 2 | = 2nNJ 2 



(Al) 



where is the number of flux quanta. We impose the 
general boundary conditions on the wavefunction for any 
particle i 



(A2) 



where r\ mn = (— i)(™ l + n +" m ) anc l we will choose 4>o = 
as our boundary condition. The operators t(a) are the 
translation operators in the presence of a magnetic fieldES 
which obey the non-commutative algebra 



t(a)t(a') = t(a + a')e i_ 5? 



(A3) 



We shall denote the two dimensional coordinate as a = 
a x + ia y and a — a x — ia y . 

In the periodic geometry, Halperin's extension of 
Laughlin's wavefunction, suitable for double layer sys- 
tems, can be written 



where 



.[Z,Z'} 



(A4) 



Vfr (mi,m 2 ,n)r / 
* re] " ' ■ ' 



[z i ,zl] = l[[ ( j>{z i -z j )}^l[[ct ) (z' i -z' j )Y 



1<J 
X 



(A5) 



Z = Y^i z i> = z' and the unprimed coordinates 
refer to electrons in the first layer while the primed refer 
to electrons in the second. This wavefunction is denoted 
(mi, m,2, n). The basic building block that we have built 
our wavefunction from is the quasiperiodic function 4>{z) 
which can be written 



(z) — w(z)exp 



(A6) 



where 
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w(z) = cxp 



6iMr) 



(A7) 



and 0i(u|t) is the odd elliptic theta function, L mn — 
K _1 (m + jiT) and r = j^-e 10 where Li = |Li||L)2|cos0. 

We have used the symmetric gauge A = |r x z in ex- 
pressing the function </>. We constrain mi and m 2 to be 
odd for fermi statistics. Note that in writing the wave- 
function in this form we have expressed the correlations 
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between the electrons but for notational simplicity sup- 
pressed the pseudo-spin part of the wavefunction which 
would properly anti-symmetrize the overall wavefunc- 
tion. This form of the double layer wavefunctions has 
been discussed previouslyll3. 

We consider a system of electrons confined to two 
parallel planes subject to periodic boundary conditions 
(PBCs) confined to the lowest Landau level. The symme- 
try analysis of this system, as introduced by HaldaneO, 
allows us to construct a Hilbert space which extracts the 
center-of-mass degeneracy as well as providing a correct 
classification of states allowing comparison with studies 
performed in other geometries. We can therefore classify 
the eigenstates of a translationally invariant Hamiltonian 
obeying [H, T(a)] = by the quantum number k defined 
to be 

r(^)l*> = (w) M exp(i^^)|*) (A8) 

where N e — Np and = Nq, allowing us to write 
v = N e /N^ — p/q. In the thermodynamic limit the states 
characterized by k = become rotationally invariant, im- 
plying that the signature of an incompressible quantum 
Hall state is a k = ground state with an energy gap to 
all excited states. As the operator ( |A8|) commutes with 
the center of mass operator T(^ sa ) = Ili^llv^) eacn 
eigenstate has a q-fold degeneracy associated with the 
action of the center of mass operator. 
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FIG. 3. Overlap of the ground state with the (3, 3, 0) state 
as a function of d/l and tunneling. Tunneling is denoted in 
the Hamiltonian as A sas and is measured in units of e 2 /47reL 
Note the change in perspective of this figure from the previous 
two figures. 

FIG. 4. Phase Diagram of v — 2/3 system as a function of 
A sas and d/l. A sas is measured in units of e 2 /4-7reZ. 

FIG. 5. Variation in energy levels as d/l is varied at 
Asas = 0. 



FIG. 6. Variation in energy levels as A sa3 is varied at 
d/l — 0. Asas is measured in units of e 2 /4nel. 

FIG. 7. Variation in energy levels as A saa is varied at 
d/l = 2.0. Asas is measured in units of e 2 /4-7reZ. 
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